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We present the results ol a search for supersymmetry with gauge-mediated breaking and Xi ~ * "/G 
in the 77+missing transverse energy final state. In 2.6±0.2 fb~ of pp collisions at y/s—1.96 TeV 
recorded by the CDF II detector we observe no candidate events, consistent with a standard model 
background expectation of 1.4±0.4 events. We set limits on the cross section at the 95% C.L. and 
place the world's best limit of 149 GeV/c 2 on the x? mass at r^o<l ns. We also exclude regions in 

the xi mass-lifetime plane for 7^0 <2 ns. 
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The standard model (SM) of elementary particles has 
been enormously successful, but is incomplete. The- 
oretical motivations [l[ and the observation of the 
'ee77+missing transverse energy (fir)' 0, candidate 
event by the CDF experiment during Run I at the Fer- 
milab Tevatron provide a compelling rationale to search 
for the production and decay of new heavy particles that 
produce events with final state photons and fir in collider 
experiments. Of particular theoretical interest are super- 
symmetry (SUSY) models with gauge-mediated SUSY- 
breaking (GMSB) These models solve the "natu- 
ralness problem" [4| and provide a low-mass dark mat- 
ter candidate that is both consistent with inflation and 
astronomical observations Since many versions of 
these models have a similar phenomenology, we consider 
a scenario in which the lightest neutralino (xi) decays al- 
most exclusively (>96%) into a photon (7) and a weakly 
interacting, stable gravitino (G). The G gives rise to 
Pt by leaving the detector without depositing any en- 
ergy [6(. In these models, the Xi is favored to have a 
lifetime on the order of a nanosecond, and the G is a 
warm dark matter candidate with a mass in the range 
0.5<to 6 <1.5 keV/c 2 Q. Other direct searches 0,0,0 
have constrained the mass of the Xi to be greater than 
100 GeV/c 2 for various points in parameter space. At 
the Tevatron sparticle production is predicted to result 
primarily into gaugino pairs, and the Xi mass (m^o) and 
lifetime (t^o) are the two most important parameters in 
determining the final states and their kinematics [l[ . Dif- 
ferent search strategies are required for xH lifetimes above 
and below about a nanosecond [ll|. 

This Letter describes a search for GMSB in which 
gaugino pairs are produced and quickly decay to the 
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77+^t+X final state, where X denotes other high-i?T fi- 
nal state particles [12j . We use a dataset corresponding to 
an integrated luminosity of 2.6±0.2 fb _1 of pp collisions 
collected with the CDF II detector ^ at ^/s=lM TeV. 
This dataset is ten times larger than the one used in our 
previous search For the first time in this channel we 
use a new photon timing system 11411 and a new model 
of the fir resolution (metmodel) [la ]. These additions 
significantly improve our rejection of backgrounds from 
instrumental and non-collision sources, which allows us 
to considerably enhance the sensitivity of the search for 
large Xi masses compared to other Tevatron searches [§] . 
We also extend the search by addressing Xi lifetimes up 
to 2 ns, which are favored for larger tti^o. 

Here we briefly describe the aspects of the detector [13 1 
relevant to this analysis. The magnetic spectrometer 
consists of tracking devices that measures the z posi- 
tion and time of the pp interaction, and the momenta of 
charged particles inside a superconducting solenoid mag- 
net. The calorimeter consists of electromagnetic (EM) 
and hadronic (HAD) compartments and is divided into a 
central part that surrounds the solenoid coil ( 1 77 1 < X . 1 ) Q 
and a pair of end-plugs that cover the region l.l<|?y|<3.6. 
The calorimeters are used to identify and measure the 
4-momenta of photons, electrons, and jets (j) [l6| and 
to provide fir information. The EM calorimeter is in- 
strumented with a timing system (EMTiming) [l4[ that 
measures the arrival time of photons. 

Our analysis begins with diphoton events passing the 
CDF three-level trigger. The combined trigger selection 
efficiency is effectively 100% if both photons have 1 77 1 < 1 . 1 
and E T >13 GeV [12j,ll5|. Offline, both photons are re- 
quired to be in the fiducial part of the calorimeter and 
to pass the standard CDF photon identification and iso- 
lation requirements [8] , with two minor modifications to 
remove instrumental and electron backgrounds ID, 17 1. 
The remaining events are dominated by SM production 
of 77; 73 witn j^lfake, and jj^jfakejfake, where -f fake 
is any object misidentified as a photon. To minimize the 
number of these events with large fir due to calorimeter 
energy mismeasurements, we remove events where the 
azimuthal angle between the fir and the second-highest 
Et photon is |A0|<O.3 or if any jet poi nts to an unin- 
strumented region of the calorimeter J15j . We require a 
primary collision vertex position with |z V crtcx|<60 cm in 
order to reduce non-collision backgrounds and to main- 
tain the projective nature of the photon reconstruction in 
the calorimeter. For events with multiple reconstructed 
vertices we recalculate the Et of both photons and fir 
values if picking a different vertex for them reduces the 
event fir- 

Non-collision backgrounds coming from cosmic rays 
and beam-related effects can produce candidates, 
and are removed from the inclusive 77 sample using a 
number of techniques. Photon candidates from cosmic 
rays are not correlated in time with collisions. There- 
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fore, events are removed if the timing of either photon, 
corrected for ave rag e path length (i 7 ), indicates a non- 



collision source [la . 117] . Photon candidates can also 



be produced by beam-related muons that originate up- 
stream of the detector (from the more intense p beam) . 
These are suppressed using standard beam halo identifi- 
cation requirements (l7j |. A total of 38,053 inclusive 77 
candidate events pass all the selection requirements. 

Backgrounds to the 77+^r final state from SM 77/ 
llfake/lfakelfake and fake fa arise due to energy mis- 
measurements in the calorimeter or to event reconstruc- 
tion pathologies. We use the metmodel [l5[ to select 
events with real and significant fa, as part of the op- 
timization, and to predict the contribution of SM back- 
grounds with fake fa due to normal energy measure- 
ment fluctuations. This algorithm considers the clus- 
tered (jets) and unclustered energy in the event and cal- 
culates the probability for fluctuations in the energy mea- 
surement to produce fa luct equivalent to or larger than 
the measured fa, P^fiuct^fa- This probability is then 

used to define a ^-significance as — log 10 (P^i ve t>fa) • 
Events with true and fake fa of the same value have, 
on average, different ^-significance. We use pseudo- 
experiments to estimate the expected ^-significance dis- 
tribution for SM events with fake fa, and the number 
of mismeasured events above a given fa -significance re- 
quirement. The jets and unclustered energy are smeared 
according to their resolution functions in the event. The 
systematic uncertainty in the metmodel is dominated 
by the uncertainty in the resolution functions. 

The metmodel does not account for reconstruction 
pathologies in SM events without intrinsic fa, such as 
a wrong choice of the primary interaction vertex or tri- 
photon events with a lost photon. To obtain the predic- 
tion for this background we model SM kinematics and 
event reconstruction using a 77 sample generated with a 
PYTHIA Monte Carlo (MC) [18| that incorporates a de- 
tector simulation [l^ . Since the pathologies from 7 j and 
jj sources are similar in nature, but not included directly 
in the simulation, we normalize the sample to the number 
of events in the inclusive 77 data sample. We subtract 
the expectations for energy mismeasurement fluctuations 
in the MC to avoid double counting. Uncertainties are 
dominated by the statistics of the MC sample, but also 
include the small differences between the measured re- 
sponse of the metmodel to MC simulation events and 
real data. 

Electroweak production of W and Z bosons which de- 
cay to leptons can also produce the 77+^r signature 
where one or more of the photons can be fake, but the 
fa is due to one or more neutrinos. To estimate the con- 
tribution from these backgrounds we use MC simulations 
normalized to their theoretical cross sections, taking into 
account all the leptonic decay modes. The Baur MC [2(| 
is used to simulate Wj and Zj production and decay 
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FIG. 1: The top plot shows the ^T-significance distribution 
for the inclusive 77 candidate sample, along with the back- 
ground predictions. The bottom plot shows the predicted Ht 
distribution after all but the final Ht requirement. 



where initial and final state radiation (ISR/FSR) pro- 
duce W/Z+jj events. The pythia MC is used to simu- 
late backgrounds where both photons are fakes: namely, 
W and Z, with photons from ISR/FSR removed, and ti 
sources. To minimize the dependence of our predictions 
on potential "MC-data" differences we scale our MC pre- 
dictions to the observed number of ej events [HI in data 
where we use the same diphoton triggers and analysis se- 
lection procedures used to select the inclusive 77 sample. 
Uncertainties are dominated by the statistics of the MC 
and ej normalization data sample. 

Non-collision backgrounds are estimated using the 
data. We identify a cosmic-enhanced sample by using 
the selected inclusive 77 sample, but requiring one of 
the photons to have t 7 >25 ns. Similarly, we create a 
beam halo-enhanced sample from events that were fil- 
tered out from our signal sample by the beam halo re- 
jection requirements [171 ] . We estimate the non-collision 
background events in the signal region using extrapo- 
lation techniques and the measured efficiencies of the 
non-collision rejection requirements [l5rj . The uncer- 
tainties on both non-collision background estimates are 
dominated by the statistical uncertainty on the num- 
ber of identified events. Figure Q] (top) shows the fa- 
significance distribution for the inclusive 77 sample, 
along with the predictions for all the backgrounds. 

We estimate the sensitivity to heavy, neutral parti- 
cles that decay to photons using the GMSB reference 
model in the mass-lifetime range, 75<m^o<150 GeV 
and T£o<2 ns. Events from all SUSY processes consid- 



6 



ered [2l| are simulated with PYTHIA followed by a de- 
tector simulation. The fraction of Xi decays that oc- 
cur in the detector volume, and thus the acceptance, de 
pend on both the lifetime and the masses of the sparti 
cles 
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The total systematic uncertainty on the accep- 
tance, after all kinematic requirements (discussed below), 
is estimated to be 7%, dominated by the uncertainty in 
the photon identification efficiency (2.5% per photon). 
Other significant contributions come from uncertainties 
on ISR/FSR (4%), jet energy measurement (2%), In- 
significance parameterizations (1%) and parton distribu- 
tion functions (PDFs, 1%). 

We determine the final kinematic selection require- 
ments by optimizing the mean expected 95% confidence 
level (C.L.) cross section limit using a no-signal assump- 
tion, before looking at the data in the signal region [22J. 
To compute the predicted cross section upper limit we 
combine the luminosity, the acceptance, and the back- 
ground estimates with their systematic uncertainties us- 
ing a Bayesian method [23]. The predicted limits are 
optimized by simultaneously varying the selection re- 
quirements for ^-significance, Ht (scalar sum of Et 
of photons, jets, and fir), and the azimuthal angle be- 
tween the two leading photons, A</>(7i,7 2 ). The large 
|5r-significance requirement eliminates most of the SM 
background with fake fr- GMSB production is domi- 
nated by heavy gaugino pairs which decay to high-i?T 
light final state particles via cascade decays. The GMSB 
signal has, on average, larger Ht compared to SM back- 
grounds so that an Ht requirement can remove these 
backgrounds effectively. Electroweak backgrounds with 
large Ht typically consist of a high-i/p photon recoil- 
ing against IF— >e^, identified as 7/ a fce|r, which means 
the gauge boson decay is highly boosted. Thus, the two 
photon candidates in the final state are mostly back- 
to-back. The SM backgrounds with fake fr and large 
Ht also have photons which are mostly back-to-back; 
the A</>(7i, 72) requirement, therefore, reduces both these 
backgrounds. 

The optimal set of requirements is slightly different 
for each point in the t^o vs. m^o space considered. We 
choose a single set of requirements to maximize the region 
where the predicted production cross section at next-to- 
leading order Q is above the expected 95% C.L. cross 
section limit. The exclusion region also takes into ac- 
count the production cross section uncertainties, which 
are dominated by the PDFs (7%) and the renormaliza- 
tion scale (3%). We find the optimal set of require- 
ments, before unblinding the signal region, to be: In- 
significance^, 7J T >200 GeV, and A^(7i, 7 2 )<tt-0.35. 
With these requirements we predict 1.4±0.4 background 
events, 0.9±0.4 of which are from electroweak sources 



(dominated by Z77 production) with real |r, 0.5±0.2 

1 ^0 001 horn non-collision 



from SM with fake fk, and 0.00L 
sources. The acceptance for 
Too«Cl ns is estimated to be 7.8±0.6%. 
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FIG. 2: The predicted and observed 95% C.L. cross section 
upper limits as a function of the x? mass at r^o^Cl ns (top) 
and as a function of the %i lifetime at m^o=140 GeV/c 2 (bot- 
tom). Indicated in green (darker shading) is the production 
cross section, along with its 8.0% uncertainty-band. In yellow 
(lighter shading) is the RMS variation on the expected cross 
section limit. 



No events in the data pass the final event selection. 
The predicted Ht distribution is shown in Fig. [1] (bot- 
tom), after all but the final Ht requirement. The data 
are consistent with the no-signal hypothesis and are well 
modeled by SM backgrounds alone. We set cross sec- 
tion limits as a function of m^o and r^o , respectively, 
as shown in Fig. O The to^o reach, based on the pre- 
dicted and observed number of events for t^o^I ns, is 
141 GeV/c 2 and 149 GeV/c 2 respectively. These limits 
significantly extend the search sensitivity beyond the re- 
sults of DO [§], expand the results to include exclusions 
for t^o<2 ns, and, when combined with the complemen- 
tary limits from CDF and LEP [lGj, |l_7| , cover the region 
shown in Fig. [3] 

In conclusion, we have performed an optimized search 
for heavy, neutral particles that decay to photons in the 
77+|t final state using 2.6±0.2 fb _1 of data. There is 
no excess of events beyond expectations. We set^ cross 
section limits using a GMSB model with Xi ~ ¥ lG, and 
find an exclusion region in the t^o-to^o plane with the 
world's best 95% C.L. lower limit on the Xi mass of 
149 GeV/c 2 at r^o<l ns. By the end of Run II, with 
an integrated luminosity of 10 fb , we estimate a mass 
reach of ~ 160 GeV/c 2 . 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
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FIG. 3: The predicted and observed exclusion region along 
with the limits found in [ToL The shaded band shows 

the parameter space where 0.5<mg<1.5 keV/c 2 , favored by 
cosmological models Q. 
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